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Abstract — Nanotransistors typically operate in far-from- 
equilibrium (FFE) conditions, that cannot be described neither 
by drift-diffusion, nor by purely ballistic models. In carbon- 
based nanotransistors, source and drain contacts are often 
characterized by the formation of Schottky Barriers (SBs), with 
strong influence on transport. Here we present a model for one- 
dimensional field-effect transistors (FETs), taking into account 
on equal footing both SB contacts and FFE transport regime. 
Intermediate transport is introduced within the Biittiker's probe 
approach to dissipative transport, in which a non-ballistic tran- 
sistor is seen as a suitable series of individually ballistic channels. 
Our model permits the study of the interplay of SBs and 
ambipolar FFE transport, and in particular of the transition 
between SB-limited and dissipation-limited transport. 

Index Terms — graphene, carbon transistors, carbon nan- 
otubes, ballistic transport, compact model, far-from equilibrium 
transport, Buttiker probes, Schottky barrier 

I. INTRODUCTION 

Since the isolation of graphene in sheets [1], [2], with their 
exceptionally promising high mobility [3], graphene-related 
materials have attracted much interest for their possible appli- 
cation in nanoelectronic devices. In particular, semiconducting 
carbon nanotubes (CNTs) [4] and single-layer or bilayer 
graphene nanoribbons (GNRs) [5] have been successfully 
employed in quasi- ID nanotransistors. 

An important issue related to carbon-based channels is the 
nature of the metallic contact at source and drain, which can 
lead to different pinning of the Fermi level and consequently 
to the formation of ohmic or Schottky contacts [6], [7]. 
The presence of SB contacts can have dramatic effects on 
device performance, because charge injection is subordinated 
to a tunneling process. However, in nanodevices with reduced 
oxide thickness, tunneling phenomena at source and drain are 
favored, and, while they often limit performance in conven- 
tional transistors, their exploitation is at the core of the concept 
of tunneling FETs [8]. 

Transport in nanotransistors is certainly far from of equilib- 
rium, but is still not fully ballistic, and currents are much 
lower than those predicted by ballistic models [9]. While 
it is perfectly clear that inelastic scattering may arise from 
the interaction of carriers with phonons and impurities, it 
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is rather complex to take into account microscopically its 
effect on transport. A powerful phenomenological attempt 
to deal with carrier relaxation and decoherence was based 
on the Buttiker virtual probes approach [10], [11], in which 
inelastic scattering is thought as localized in special points, 
spaced by a length defined as "mean free path". The Buttiker 
approach was also introduced in microscopical models based 
on tight-binding Hamiltonians [12], and recently extended to 
deal, via a quantum Langevin approach, with ID conduc- 
tors [13]. In [14] the Buttiker probes approach to inelastic 
scattering was employed in a simulation, based on the non- 
equilibrium Green's functions formalism, of a non-ballistic 
silicon nanowire transistor. 

Fully microscopical analysis of inelastic scattering due 
to specific mechanisms such as phonon scattering, with the 
non-equilibrium Green's functions approach, has also been 
addressed by adding a proper self energy correction on a site- 
representation propagating Hamiltonian by Jin et al. [15] and 
by M. Gilbert et al. [16], [17]. 

As far as analytical models are concerned, transport in 
quasi- ID FETs is generally treated as purely ballistic or with a 
drift-diffusion assumption as in Ref. [18], [19], [20]. A largely 
invoked approach to treat partially ballistic transport including 
the effects of backscattering was proposed by Lundstrom et 
al. [21]. This approach, that is easily included as a correction 
to ballistic models, has the merit of offering a very simple 
and synthetic picture but does not allow a full description 
of the seamless transition from ballistic to quasi-equilibrium 
drift-diffusion transport. Recently a rigorous semi-analytical 
model based on the Buttiker virtual probes approach [10], [11] 
has emerged, in which a non-ballistic transistor is seen as a 
suitable chain of N ballistic channels, where N is the ratio of 
the channel length to the mean free path, or equivalently as a 
series of drift-diffusion and a ballistic FET [22], [23], [24]. 

In this work we propose a semi-analytical model based on 
the virtual probes approach, which describes one-dimensional 
FETs, treating on equal footing Schottky barrier contacts and 
FFE transport conditions. In Section II we summarize the gen- 
eral analytical description of graphene nanoribbons subbands, 
density of states, equilibrium charge density, extensible also 
to the carbon nanotube case. In Section III we present a WKB 
approximation of the tunneling probability through Schottky 
barrier contacts, yielding analytical expressions for the trans- 
mission based on two different levels of approximation for the 
energy dispersion curves of GNRs (or CNTs). In Section IV 
a model for a single ballistic transistor with SB contacts is 
presented, compared with data from numerical simulations. In 
Section V we propose a compact model, based on Buttiker 











sio 2 






S 




D 






Si0 2 















Fig. 1 . Sketch of a DG-GNR field effect transistor, considered as an exemple 
for the application of our model. 

virtual probes approach able to deal with both intermediate 
transport and SB contacts, and use it to study the interplay of 
SB and dissipative transport. 
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Fig. 2. Comparison of the subbands of a A-GNR with 12 dimer lines 
between a numerical tight-binding calculation and our analytical result with 
edge corrections. Valence bands are symmetrical. 



II. DISPERSION RELATION AND DENSITY OF STATES 

The dispersion curve of an armchair GNR with N dimer lines 
can be obtained analytically by cutting techniques, analogous 
to that used for CNTs in [25], from the 2D graphene tight- 
binding dispersion. The subband dispersion curves correspond 
to ID segments of the graphene Brillouin zone with the con- 
fined wavevector quantized as k a = with a - 1,2, . . .,N. 
The dispersion curve of the subband a, referred to midgap, is 



E a {k) = ±V\\ 



4 cos — - — A„ + 4A" 



1/2 



(1) 



energy dispersion relation, which we will refer to as the full 
band (FB) approximation when applied to FET modelling, is 
therefore 

E c a (k) = E a (k) + 8E a (k). (4) 

The comparison between numerical tight-binding calculations, 
with edge effects taken in account, and the analytical result 
with perturbative corrections, for a A-GNR of 12 dimer lines 
is shown in Figf2] The agreement is very good, especially at 
k = 0, where © reproduces the results of [27]. For simplicity 
we define here the band edges as s a = £^(0). 



with A a = cos(-^j-). We note here that a dispersion relation 
totally analogous to (Q]) applies to zig-zag (N,Q) CNT, with 
the only difference that in the place of A a we have to use 
A™ T = cos (f ), where a is the subband index of CNTs [26]. 

Therefore, much of the results for GNRs obtained here and 
in the following of the paper, with the exclusion of the edge 
corrections, can be directly generalized to the zig-zag CNT 
case. The edge of the a-th subband E a (0) is expressed as 



E a {U) = ±V{\+2A a ). 



(2) 



Let us note that the lowest lying subband is given by the 
value of a for which A a + j is minimum. The edges of the 
nanoribbon are laterally exposed to vacuum and experience a 
different chemical environment, therefore the hopping param- 
eter between carbon atoms at the edges tends to be slightly 
different. We can, at least partially, account for the presence of 
edges via a perturbative approach to the first order [27]. The 
perturbation theory to the first order leads to the following 
eigenenergy corrections: 

SE a (k) = (±) a H e a d a = (±)«^j sin2 (^^) cos (^)> (3) 

with v = 0.12 eV the energy correction of the hopping 
parameter at the edges in the tight binding Hamiltonian. The 
correction has a positive or negative contribution depending on 
the wavefunction parity with respect to the two asymmetric 
carbon atoms, which are connected by the edges. Therefore 
if A a > —j we have a positive contribution (±) a = 1, 
otherwise a negative one (+) a = -1. The edge-corrected 



A. Approximated expressions 

In modelling nanotransistors only the lowest laying subbands 
matter, in which the relevant transport phenomena take place. 
For these lowest lying subbands often an effective mass (EM) 
approximation is invoked 

e„ + ~z~r~r~ 5 (5) 
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in this case the following effective mass for the a-th mode can 
been employed 

2 h 2 s a 



Ma ~ 3a^V 2 A a - 
The DOS in EM approximation is given by 
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with E expressing the 'kinetic energy', i.e. the energy calcu- 
lated with respect to the band edge s a . 

The EM approximation is rather crude, and an intermediate 
(I) approximation, between the FB and the EM, can be the 
use of the dispersion curve 



E' a (k) = 
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Fig. 3. Energy dispersion curve, and the corresponding density of states, of 
the lowest conduction subband of a A-GNR with 12 dimer lines. A numerical 
tight binding result is compared with our FB analytical result and with EM and 
I approximate dispersions. In the energy range considered here the agreement 
between numerical, FB and I approximations is excellent. 



In Fig|3]we compare the lowest band dispersion curve and the 
corresponding DOS for a GNR with 12 dimer lines. Both the 
FB and the I approximations reproduce quite well numerical 
tight-binding calculations, and give similar DOS, of course 
the I dispersion is only accurate for energies £«V. The EM 
approximation instead remains quite accurate only for about 
£<0.1eV. 



B. Carrier density 

The carrier density affects both electrostatics and transport 
properties. Here we develop a similar analysis to what done 
in [26] for carbon nanotubes. The electron carrier density per 
subband can be expressed as 



I 



/( 



E + e a - q<pc - n 
kf 



)D a (E)dE, 



(10) 



where fi is the Fermi level, f(x) = (1 + expjc) -1 the Fermi- 
Dirac distribution, and e'° p is top edge of the a-th subband, 
that for most purposes can be taken as oo, due to the finite 
extension of f(E). (p c is the electrostatic potential in the device, 
which rigidly shift the levels. Because the non-negligible 
contribution to dTOb comes from states near s a , we can use the 
intermediate expression for the DOS D[ y (E). If we consider 
a non-degenerate situation (s a - 3kT > fi), typical of sub- 
threshold regimes in FETs, we obtain 
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with z = E/s a and B = (kTy 1 . With a partial integration and 
recognizing the modified Bessel function of the second kind 
K\, the charge density can be expressed as 



n a = — y/M a s a 



(12) 



In order to give an estimation of the Bessel function K\ which 
has no closed form, we can adopt the approximation [26] 



Ki(x) 



K 1/2 (x) + K V2 (x) 
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arriving in the end to express the charge density as 



M a 1 + 2Bs a 



(13) 

(14) 
(15) 
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with essentially the same form of 3D bulk semiconductors. 

III. TUNNELING OF SCHOTTKY BARRIERS 

Our aim is to provide an analytical description of the tunneling 
through SB contacts. The first step is to model in the simplest 
way the potential decay occurring near the source and drain 
contacts. The potential inside a transistor channel is described 
by the a 3D Poisson equation 



V 4>(r) = 

e 



(16) 



together with the boundary conditions enforced by voltages V s , 
Vd, V g at the source, drain and gate leads. In the evanescent 
mode analysis approach the electrostatic potential inside a 
nanotransistor <p(r) is thought as the sum of a long-channel 
solution 4>i{r), which satisfies the vertical electrostatics, plus a 
short-channel solution <fi*(r), called evanescent mode, responsi- 
ble of the potential variation along the channel [28]. The short- 
channel solution is obtained, solving the Laplace equation for 
the device with an adequate expansion in harmonic functions. 
As a matter of fact the short-channel solution near the source 
contact results in an exponential profile 

f{r) cc R(^)e- Z/A . (17) 

where R(r\) describes the solution in the channel cross section 
and A comes to be a the natural scale length for the potential 
variation in the device. The actual value of A depends on 
the details of the device geometry, however in double-gate 
(DG) configuration, and considering that in general, in carbon- 
based FET, the oxide thickness is significantly larger than the 
channel thickness, the asymptotic value A = (2t ox + t c h)jn can 
be assumed. In the case of a cylindrical GAA-CNT FET, an 
explicit calculation of A via evanescent mode analysis has been 
performed in [29]. 

We follow this line and assume that the channel potential 
rigidly shifts the confinement eigenvalues s a , where a runs 
on the different subbands. Now we are interested only in the 
potential inside the restricted zone of the graphene channel 
<pc(z), in which it can be assumed as a constant (which is 
strictly true in subtreshold regimes), and we consider its 
variation only along the channel direction. The long channel 
solution inside the channel is reduced to <pL(r) ~ <pc, where <p c 
is solely imposed by the vertical electrostatics, while the short- 
channel solution has the form ST% . Therefore the potential in 
the channel <f> c (z) can be expressed as 



c (z) = 4>c + —e~ z,A 

q 



(18) 
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with <p c - <p(oo) fixed by the vertical electrostatics and A s 

imposed by the boundary condition at the SB contact A s = 

is) 

Egg — Sl + q<pc where L refers to the lowest lying subband, 
due to the Fermi level pinning at the metal/semiconductor 
interface. E SB is the Schottky barrier height on the first 
conduction subband with respect to the source Fermi level. 
The charge injected from the source with energy lower than the 
barrier have to tunnel in order to reach the channel. We need 
to calculate the transmission through an exponential decaying 
barrier of the kind 



t=1.5nm 



t = 4 nm 



E SB (z) = A s e-*' A 



(19) 



with the height A s dependent on the electrostatic potential <p c . 
We note that however that if the band bending exceeds the 
energy gap 2s a , carrier with energy < E < A s — 2e a , will 
experience a SB of an height A s — E + 2e a . 

In order to estimate the behavior of a nanotransistor it is 
essential to accurately describe tunneling phenomena, both in 
traditional FETs and in TFETs. In this section we compare the 
tunneling calculated with WKB approximation in a full band 
approach (FB-WKB), within the effective mass approach (EM- 
WKB) and intermediate approximation (I- WKB). FB-WKB is 
more complex to implement and requires a numerical solution 
of the integral 



ln(T(E)) 
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I[k z (z; E)]dz 
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While for the others two an analytical expression for the 
tunneling T(E) can be obtained. 

A. Effective-mass WKB approximation 

The transmission coefficient obtained via WKB approxima- 
tion is given as 



T(E)= e x ' £<A < (21) 

1; E>A S 

where z\,zi are the classical turning points are 

\ E 

z\ = 0; zi - -A In — 
LA, 

The transmission coefficient can be analytically calculated in 
In T(E) = -4/l-i/ 1 - J - tan" 



ti 2 



A s -E 



(22) 



B. I WKB approximation 

Let us consider a dispersion curve of the kind ([8j. The 
turning points with a barrier like (fT9l are the same as d22l) . but 
now, under the barrier, the imaginary part of the wavevector 
as a function of the energy is given by 



I[k z ,E] 
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with 
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Fig. 4. Transmission probability for an electron in the lowest laying subbands 
of a A-GNR with 12 dimer lines. The FB result comes from the numerical 
integration of (28) , while we have analytical transmission probability in I (24} 
and EM 1221 approximation. 



The integration ( l20l i. for E < A s , leads to the WKB tunneling 
probability 
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where we introduced the abbreviations 



Ri = y/a 2 - (b- l) 2 ; R 2 = ib 2 - a 2 . 

C. Full-band WKB approximation 

For an armchair GNR, subband dispersion curves are in the 
form (U), from which we can express the wavevector as a 
function of the energy as 

2 



fV3 



: arccos x, 



(25) 



with the substitution u = -z/A and normalizing all quantities 
to A, x given by 



(f +5-e u ) 2 -a 



(26) 



where we introduced 



6 = 



a - 



(l+4A ff )V 2 



4A„V 2 
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In the integration domain of ( 1201 ). the argument x of the inverse 
cosine function has module larger than 1, and therefore 

2 



I(k) 



<V3 



ln z + yjz 2 + 1 



leading to the WKB tunneling probability 

ln£/A s 



InTXE) = — 
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In be + Vx 2 + 1 \du. 



(27) 



(28) 



In Figj4] we compare the tunneling coefficients, calculated 
with the EM, I and FB WKB approaches, for a SBs of 
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height 0.5 and 1 eV, and for a A typical of DG A-GNR with 
t ox = 1.5 and 4 nm. Essentially the intermediate approxi- 
mation completely reproduces the FB tunneling probability, 
while a significant deviation is observed with the EM-WKB 
approximation for E < 0.5A S . Therefore the intermediate 
approximation seems an optimal approximation for compact 
models in order to reduce the computational times retaining 
high accuracy. 

IV. SCHOTTKY BARRIER BALLISTIC FET 

We consider here a ballistic transistor with Schottky barrier 
contacts at source and drain, as shown in Fig|5] As usual 
in compact models, we assume a complete phase random- 
ization along the channel, neglecting phase resonances in the 
transmission probability of the two tunneling barriers, while 
multiple reflection events are taken into account. Between 
two tunneling barriers, the forward and backward distribution 
functions are modified by the multiple elastic scattering [30], 
[29]. The overall mobile charge, given by the sum of forward 
and backward going charge carriers in the channel, can be 
expressed as 



a & 1 T C SB thermionic 



T d (2-T s ) ,. 

.)+ — - — fii a4 ) 



with i - e, h for the electron and hole charge, where 
_ E - q<p c +s a - fi sW 

h _ Hs(d) - E + qcf> c - s a 

n a ,m - hT 

r = T, + T d - T s T d , 



and 



(29) 

(30) 
(31) 

(32) 



where T s , T d are the tunneling coefficients at source and drain, 
depending on both energy and channel potential. In order to 
compute the channel potential <f> c , and, through it, the subband 
energies, the total mobile charge Q - Qh — Q e must be equal 
to the charge induced by the electrostatic coupling of channel 
with gate, source and drain through the capacitances C g , C s , C d 
respectively: 

Q(0 C ) = - Yj CiiVi-Vfsj-fa), (33) 

i=g,s,d, 

where Vfb j = 4>i ~ Xg is me flat band voltage, given by the 
difference beween the contact workfunction and the graphene 
electron affinity. 

The current is obtained with the Landauer-Biittiker formal- 
ism, which, accounting for the tunneling, takes the following 
form 



h(4>c) 



c a —t 



T s T d 



U(i a J-f(j a4 )]dE. (34) 



with i - e, h accounting for the current of electrons and holes, 
and the total current given by / = I e — //,. We note that (|29l 
and d34l i include both tunneling and thermionic contributions. 
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Fig. 5. Conduction band edge profile of a SB nanoscale FET. The thermionic 
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Fig. 6. Comparison between our model and numerical simulations from 
ref.[31]. The transfer characteristics (a) at V ds = 0.5 V and characteristics 
(b) at V g = 0.75 V of a ballistic double-gate armchair GNR FET, with Ohmic 
and SB contacts of height E g /2 as 0.3 eV are shown. Assuming a GNR 
thickness of about 1 nm we obtain A as 1.3 nm. 



We apply our model to the case of a double gate armchair 
graphene nanoribbon transistor (DG A-GNR FET) with both 
Ohmic and Schottky barrier contacts. In Fig|6] we compare 
the transfer-characteristics (a) and the output characteristics 
(b) of a ballistic armchair GNR FET, obtained with our 
model and with numerical simulations based on the non- 
equilibrium Green's function formalism in Ref. [31]. The S1O2 
gate oxide thickness is 1.5 nm, the armchair GNR lattice is 
characterized by 12 dimer lines, which correspond to a width 
of 1.35 nm and a bandgap of 0.6 eV. We employed here 
the intermediate analytical description of the GNR subbands 
and density of states ((9]). The source and drain capacitances 
C s , C d are introduced because of the short-channel nature of 
the GNR simulated in [31] and are fixed, with respect to the 
gate capacitance C g = 1.1 X 10~ 10 F/m, to C s = C d = 0.1C g . 
The agreement between the numerical simulations and our 
compact model, for both curves (Figj6j> with ohmic and SB 
contacts, is very good, demonstrating that the effects of SBs 
are well accounted for. 

V. SB TRANSISTORS IN INTERMEDIATE TRANSPORT 
REGIME 

To describe dissipative transport, we follow the approach 
developed in [22], [23] for a 2D MOSFET for the non- 
degenerate and degenerate cases, and in [24] for quasi- ID 
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Fig. 7. Chain of N ballistic transistors with SB contacts at source and drain 
(first and last contact). As explained in the text, the chain of ballistic transistors 
can be described as the series made by a central DD section accounting for 
dissipative transport in the N—2 internal nodes and by head and tail ballistic 
transistors accounting for the SB contacts with source and drain. 



FETs. Such treatment is here expanded to include ambipolar 
devices. We recall that, within the Biittiker probes approach, 
inelastic scattering is thought as localized in special points, 
spaced by a length defined as "mean-free path" i. The virtual 
probes act as localized reservoirs along the channel, in which 
carriers are fully thermalized in equilibrium with the probe 
quasi-Fermi energy fi„. Transport from one virtual probe to the 
next is considered purely ballistic. We have a drift-diffusion 
transistor when the channel length is much longer than the 
free mean path, that from our point of view it is equivalent to 
have a long enough chain of ballistic transistors, as rigorously 
shown in [22]. On the contrary, when the number of internal 
contacts is small, transport is far-from-equilibrium, and fully 
ballistic in the limit N = 1. 

A transistor with SB contacts in the FFE transport regime is 
therefore modeled as a series of individually ballistic channels, 
connected by fully thermalizing virtual probes placed at x„ = 
nl with n = 1, . . ,,N — 1, with electrochemical potential 
Head and tail of the series are connected to source and drain 
through SB contacts as sketched in Fig|7J and boundaries are 
fixed as jjQ = /u s and = fid- In the n-th ballistic channel yu„_i 
and /j„ act as source and drain, simultaneously solving (|29| | 
and d33l we can fix the channel potential $ • In me same 
manner the current in the n-th channel is obtained with d34t , 
imposing /i„_i and fi„ as source and drain Fermi levels. Since 
the current /„ in any n = l,...,N FET must be equal to Ids, 
we have N equations determining the local Fermi energies fj. n . 
We note that a distinction between ballistic internal channels 
(B) and boundary channels with source (B s ) and drain (Bd) 
can be made. In fact the first and the last ballistic channels 
are characterized by SB contacts with metallic source and 
drain, while internal channels, in the region between the 
fictitious virtual probes, can be treated as ohmic transistors. 
The numerical solution of the complete chain of N elements: 
2 of the boundary kind and N - 2 of internal kind, will be 
addressed as the B(N) model. 

Now we note that for the internal part of the chain the 
analysis developed in [24] applies. In particular it has been 
shown that the current in an ohmic-contact ballistic chain of 



N elements, after a linearization procedure, can be arranged 
in a to a drift-diffusion-like form (that we refer as the DD(N) 
model) in which the current is calculated through the formula 

v " 

q 2 r(l){ 



Ids = 



uhL 'Z \{ F -MV\)-F-i(il h a [V]))dV, 



(35) 



where F-i(x) is the Fermi-Dirac integral of order -1, T the 
gamma function and 



= (q<p c - qV - s a )/kT, 
(-q<f> c +qV + S a )/kT. 



(36) 
(37) 



We note that 77, not only directly depends on V, but also 
indirectly through <p c , which is self-consistently imposed by 
the linearized vertical electrostatics 

Qm = C g (V g - V FB - <P C [V} ) 

Qm = X a yp§^(F 1/2 (if a [VT) - F 1/2 (77*[V])}. 

The linearized DD model d3~5T l has also the advantage of 
dealing with non integer = L/(, and is therefore more 
flexible than the ballistic chain itself. As noted in [24], (l35l ) 
can be rearranged in a local form, analogous to a DD equation 
la = HaQa jf, where the degenerate mobility (we consider now 
a monopolar regime) is given by 

qv a { F-i[tf a ] 



2kT F. m [rf a ] 



(38) 



with v a - yj^- the mean carrier velocity. This expression 
gives us a link between N - L/( and the mobility. 

We can now model a SB transistor in intermediate transport 
regime as a series of B S -DD(N)-B S segments, with two nodes 
between the boundary channels and the internal segment, 
characterized by electrochemical potentials that can be fixed 
exploiting the current continuity in the device. We will refer 
to this macro-model as the BDDB(N) model. This compact 
model permits to analyze both the presence of Schottky barrier 
contacts and far-from equilibrium transport condition, while 
keeping low the computational burden, especially with respect 
to numerical simulations including dissipation. 

We now analyze the effects of inelastic scattering on the 
performance of a DG A-GNR FET. In non-ballistic transport 
(increasing AO the transfer characteristics (Fig|8]l vertically 
shift, in a semilog plot, as expected due to the mobility 
reduction. It is interesting to note that the effect is more marked 
in the subthreshold region and, consequently, an increase of 
the I on /I ff ratio as a function of Af is observed, as shown in 
the inset. In ballistic models with positive Vds, in subthreshold 
regime, tunneling from the drain leads to hole accumulation 
under the channel, which increases the quantum capacitance 
and reduces the control over channel. Subsequently a larger 
subthreshold swing and a lower I on /I ff is obtained. 

An accurate analysis of the SB effects on output characteris- 
tics can be performed calculating the differential conductance 
g = dlds/dVds- In fig |9] we compare the output characteristics 
and the differential conductance for a device with t ox = 5 nm, 
with a SB height SB= 0, 0.25, 0.5 E g . Note that the presence 
of Schottky barrier contacts is more relevant in transistors with 
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Fig. 8. Transfer characteristics of a ballistic chain made of a series of N DG 
GNR FETs, with N= 1,2,5 and V rfj = 0.5 V, calculated with our model. In 
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Fig. 9. Output characteristics and differential conductance for of a GNR 
ballistic chain made of N DG GNR FETs, with N = 1, 2, 5, 10, 20. Three 
different pinning of the SB with respect to the conduction band are considered: 
without Schottky barrier (SB=0), with a SB of E g /4 and E g /2. The details of 
the device are the same as Fig[6] except t ox = 5 nm, for which all the features 
due to the SB are enhanced due to the thicker barrier. 



a looser vertical confinement, where the tunneling barriers 
are thicker. We observe that in samples with SB= eV 
the output characteristics concavity is always negative, and 
the differential conductance is monotonously decreasing with 
Vds- If the SB height is finite the differential conductance 
acquires a non-monotonous behavior, which well describes 
the "S" shaped concavity change of the characteristics curves 
before reaching saturation, especially evident in thicker SB 
devices. It is interesting to note as, apart from a reduction 
of the maximum saturation current, larger ballistic chains 
(larger AO, in which a higher inelastic scattering is active, 
lead to a smoothening of the non-monotonous dependence of 
g on Vds- In this fact we can recognize a gradual transition 
between devices in which the characteristics are dominated 
by SB contacts and devices in which inelastic relaxation is 
predominant. 

In electron-hole symmetrical materials as undoped graphene 
nanoribbons or carbon nanotubes, the relative SB height with 
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Fig. 10. Transfer characteristics of GNR devices increasing t ox = 1, 3, 5 nm, 
calculated for Vj s = 0-5 V. Ballistic chains of N = 1, 2, 5, 20 are drawn, for 
devices without Schottky barrier (SB=0), with a SB of E g /4 and E g /2 are 
shown. Arrows indicating the shift of the transfer characteristics curves with 
TV are also added as a guide for the eyes. 



respect to the bandgap determines the position of the minimum 
of transfer characteristics, it influences their shape and their 
symmetry (see FigJTOt. A SB of height E g /2 preserves the 
bandstructure electron-hole symmetry and therefore results in 
transfer characteristics which span symmetrically from the 
current minimum off state (placed at V g = Vj s /2). Curves 
calculated reducing the SB height for electrons (for E g /4 and 
0) show a growing asymmetry, with weaker hole currents 
and larger electron currents, together with a shift of the 
transfer characteristic minimum to lower values of Vg. This 
phenomenon is prominent in thicker SB devices such as the 
t ox = 5 nm FET, but well observable also in a t ox = 1 nm 
device. The increase of the lateral confinement leads in fact to 
an almost linear increase of the SB thickness and therefore 
all tunneling processes become harder. As expected, if we 
increase the dissipative phenomena (increasing N) a reduction 
of the current is observed. But more interesting, while the 
SB= Eg/2 curves vertically shift along the segmented line, 
the shift of the other curves is diagonal, note in fact the 
horizontal shift of their minima with N. Moreover, increasing 
N, the minima seem to converge towards the value V g = Vd s /2, 
typical of a symmetrical ambipolar device. This is yet another 
signature of the growing importance of inelastic transport over 
the SB contacts. Therefore, for sufficiently well-confined FET, 
we can expect in quasi-ballistic GNT/CNT devices to clearly 
observe a SB behavior, which become more and more subtle in 
dissipative regimes. To quantify the relative importance of the 
Schottky barrier in determining the symmetry of the transfer 
characteristics we made the following physical estimation: 
SB= 0.5 E g corresponds to the symmetrical case, therefore if 
we impose a different SB the change in the conductance will 
be exponential in the SB difference £>£sb as 



-'sb - e\p-(-^(2m<5£ , sB) 1/2 



(39) 



as can be obtained estimating the differential conductance of 
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a device with a SB source contact at the source Fermi level. 
This quantity is in fact dominated by the tunneling coefficient 
( l22l >. This difference in the conductance is relevant as long as 
it is greater than the conductance due to the DD(N) chain. We 
obtain 



0.1 V 



T = — * Nexpl-t,, * ^-(2m6E SB ) l/2 
8n [ nn 



(40) 



Employing this formula we can calculate the N = N s corre- 
sponding to y = 1 for different SB value and oxide thickness, 
as shown in the following table 



N s 


tax = 1 




t ox = 5 (nm) 


SB= 





4 


70 


10 a 


SB= 


0.25£ g 


7 


450 


10 4 



N s gives a rough estimation to the number of nodes (i.e. L/( 
ratio) needed to make the transfer characteristics symmetrical, 
in spite of the presence of a SB. As can be observed comparing 
these values with the behavior of curves in Fig|T0] the t ox = 
1 nm curves with SB= 5 and 10, respectively for SB= Q.25E g 
and 0, are quite symmetrical in accordance with N s = 4 and 
N s = 7 found by our calculation. The minimum of the curve 
N — 10 with SB= 0.25£g comes near to the symmetrical 
values, but still misses it being our estimation N s = 70. Other 
curves are highly asymmetric being N <s N s . 

A typical parameter used to characterize the transport 
regime in quasi-ballistic devices is the ballisticity index 
Bindex - I/h, which is the ratio of the actual current to 
the current corresponding to an analogous device in a purely 
ballistic transport regime (N = 1). In FigQT|we analyze the 
role of the SB contacts in determining the ballisticity index as 
a function of N, and therefore as a function of the degree of 
inelastic relaxation. In general to lower SB heights correspond 
a faster variation of the B in( j ex with N, with a sudden drop of 
the ballisticity as function of the number of nodes, after which 
a slower decrease is observed. SBs affect in particular the 
ballisticity index calculated for lower Va s , due to the concavity 
of the output characteristics, while larger source-drain voltages 
reduce the relative importance of SB with respect to inelastic 
mechanisms. Calculations with t ox = 5 nm reveal the increased 
importance of SB contacts, and reflect the presence of the 
inflection in the output characteristics, with a concavity change 
before saturation. In particular, for higher value of the SB we 
observe a slower dependence of the B int j ex on N, because the 
current is calculated in a bias point of the characteristic curve 
of strong "s" curvature. Physically, it means that the current 
flowing in the device is mostly limited by the injection through 
the tunneling barriers. 

VI. CONCLUSIONS 

We have presented a semi-analytical model dealing with 
ambipolar one-dimensional Schottky barrier transistors in in- 
termediate transport regimes between fully ballistic and quasi 
equilibrium, i.e. governed by the drift-diffusion model. We 
have introduced simplified, but accurate, descriptions of the 
Schottky barrier profiles and of the electrostatics, and analyti- 
cal approximations of the tunneling coefficients of the Schot- 
tky barriers. We demonstrate that a Schottky barrier transistor 
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Fig. 11. Ballistic index of a ballistic chain of N elementary GNR FETs with 
SB is varied from to Eg/2 as a function of N. The source-drain voltage is 
set to V ds = 0.1 V and 0.5 V for V, = 0.75 V. 



can be modeled as three transistor is series, with common 
gate voltage. The central one is a drift-diffusion transistor, 
with mobility dependent on the degree of degeneracy of the 
one-dimensional carrier gas. The other two transistors are 
ballistic FETs with a Schottky barrier contact corresponding 
to the external actual contacts (source or drain). In the case 
of ballistic transport, our model allows us to reproduce the 
results of a 3D numerical Poisson-Schroedinger simulator. In 
the case of very long channel, with respect to the mean free 
path, current is limited by the central drift-diffusion transistor. 
The model allows very directly to investigate the transition 
from barrier-limited transport to channel-limited transport. 
Our semi-analytical model represents an accurate and simple 
way to gain physical insights into the behavior of nanoscale 
transistors with Schottky barrier contacts, including most the 
relevant physics at a very low computational cost. 
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